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CHAPTER  1 


INTRODUCTION 


This  report  discusses  the  setup  and  operation  of  a  .computerized 
apparatus,  shown  In  Figure  1,  for  use  In  measuring  dielectric  properties  of 
selected  materials  under  Investigation.  These  measurements  are  directly 
applicable  to  microwave  Interferometerlc  measurements  where  the  microwave 
wavelength  In  the  material  under  investigation  is  a  variable  dependent  on  its 
dielectric  constant,  the  microwave  test  frequency  and  the  boundary  conditions 
of  the  material.  The  procedure  discussed  here  Is  based  on  the  use  of  a 
slotted  line  to  measure  the  parameters  necessary  for  the  computation  of  the 
dielectric  properties  whose  boundary  conditions  are  a  rectangular  waveguide. 
The  accuracy  of  the  dielectric  constant  computation  is  to  a  large  extent 
dependent  on  the  measurement  of  the  location  of  a  null,  or  voltage  standing 
wave  ratio  (VSWR)  minimum,  of  the  interference  pattern  resultant  from  the 
•  incident  signal  and  the  reflected  wave  from  the  sample.  A  stepper  motor  with 
movement  accuracy  of  ±1  micrometer  was.  applied  to  improve  the  precision  of  the 
spatial  measurements  compared  to  the  manual  technique  previously  used  to 
locate  the  voltage  minima.  This  use  of  the  stepper  motor  also  improved  the 
VSWR  dramatically  and,  coupled  with  the  Improvement  in  minima  location, 
enhanced  the  accuracy  of  both  the  real  and  imaginary  parts  of  the  dielectric 
constant  computation. 

In  November  of  1972,  the  Agricultural  Research  Service  of  the  U.  S. 
Department  of  Agriculture  published  a  report1  outlining  a  computer  program 
that  had  been  developed  for  calculating  dielectric  properties.  This  program, 
originally  written  in  FORTRAN,  had  been  adapted  to  run  in  HP-BASIC  on  an 
HP  9845B  computer.  With  the  computer  program  and  a  precision  stepper  motor 
available,  it  was  possible  to  incorporate  these  resources  into  a  single 
apparatus  capable  of  automated  dielectric  property  measurements.  The  initial 
objective  of  these  automated  measurements  was  to  facilitate  research  into  the 
development  of  a  material  with  equal  permittivity  (dielectric  constant)  and 
permeability  for  dielectric  antennas.  However,  the  procedures  that  have  been 
developed  should  be  suitable  for  use  wherever  dielectric  constant  measurements 
are  needed. 
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CHAPTER  2 


EXPERIMENTAL  ARRANGEMENT 

The  dielectric-properties  computerized  setup  is  based  on  the  use  of  a 
slotted-line  probe  carriage  and  a  two-point  method  for  measuring  the 
dielectric  constant  involving  the  solution  of  a  transcendental  equation.2,3-*  A 
tuned  probe  mounted  in  the  carriage  is  connected  to  an  HP-415D  SWR  meter.  The 
carriage  allows  for  the  positioning  of  the  probe  to  locate  the  standing  wave 
minima  while  the  probe  slides  along  a  relative  scale.  The  carriage  scale 
gives  position  measurements  which  are  used  in  the  calculation  of  the  microwave 
parameters  and  allows  measurements  to  be  taken  to  within  +.1  mm  over  a  range 
of  130  run.  This  measurement  can  be  Improved  by  mounting  a  dial  gage  onto  the 
carriage  which  enables  measurements  to  be  taken  to  within  +.01  mm,  but  only 
over  a  range  of  approximately  25  mm.  This  small  range  is  inconvenient  for  the 
determination  of  complete  dielectric  properties.  In  place  of  this,  a  stepper 
motor  was  connected  to  the  probe  carriage  which  would  drive  the  probe  over  a 
120  mm  range  with  an  accuracy  of  ±.001  mm.  The  addition  of  the  stepper  motor, 
with  a  controller  allowing  a  computer  to  regulate  its  movements  via  an  IEEE- 
488  interface  bus,  made  it  possible  to  automate  the  measurement  acquisition 
process. 

An  HP-8350A  sweep  oscillator  with  an  IEEE-488  interface  was  available  as 
the  microwave  source.  This  source  provides  the  option  of  square  wave 
modulation  of  the  signal,  which  is  necessary  to  obtain  gain  from  the  narrow 
band  amplifier  of  the  standing  wave  ratio  (SWR)  meter  used.  Other  necessary 
features  of  this  source  are  frequency  control  and  power  level  control  via  the 
IEEE-488  interface.  The  frequency  range  with  the  available  sweep  oscillator 
plug-in  was  .01  to  8.4  GHz  initially,  so  when  working  in  X-band  (8.2-12.4  GHz) 
waveguide  as  we  did,  the  usable  frequency  range  was  then  only  8.2  to  8.4 
GHz.  A  plug-in  that  will  totally  cover  the  X-band  frequency  range  has  been 
employed  only  recently  and  was  used  in  some  measurements  that  are  reported 
here. 


An  HP-3438A  digital  multimeter  with  IEEE-488  interface  was  connected  to 
the  voltage  output  of  the  SWR  meter.  A  voltage  reading  is  sent  to  the 
computer  corresponding  to  the  power  level  measured  by  the  SWR  meter.  A  full- 
scale  reading  outputs  500  mV  to  the  voltmeter  with  0  mv  output  for  a  minimum 
scale  reading.  This  remains  constant  for  any  scale  used  on  the  SWR  meter, 
although  the  SWR  meter  remained  only  on  its  most  sensitive  scale.  A  voltage 
reading  received  from  the  SWR  meter  on  the  multimeter  is  output  to  the 
computer  and  then  checked  before  the  power  of  the  source  is  adjusted.  This 
way,  the  power  seen  by  the  tuned  probe  can  be  held  low  with  the  voltage 
reading  maintained  below  100  mV,  thus  keeping  the  probe's  crystal  operating  in 
its  square-law  (linear)  region.  This  also  avoids  a  temperature  change  or 
humidity  bake-out  due  to  high  power  dissipation  in  the  crystal  detector.  When 
the  probe  is  located  at  a  standing  wave  minimum,  for  example,  a  reading  of  100 


mV  will  correspond  to  the  power  level  of  the  source  being  raised  to  its 
highest  level.  The  concept  of  having  maximum  power  output  from  the  source 
when  the  probe  is  located  at  a  standing  wave  minimum  is  a  concept  that  is  very 
important  for  the  understanding  of  the  automated  program  shown  in  Appendix 
A.  The  principle  of  operating  in  the  square  law  region  of  the  probe's  crystal 
is  a  basis  upon  which  this  program  is  developed. 
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CHAPTER  3 

PROGRAM  DISCUSSION 

The  program  is  begun  with  the  empty  sample  holder  mounted  in  place  on  the 
microwave  circuit.  This  is  done  to  obtain  a  base-line  measurement  for  the 
circuit  with  the  sample  holder  to  be  used.  The  program  asks  for  the  initial 
room  conditions  of  temperature  and  relative  humidity  as  well  as  the  frequency 
to  be  used.  The  sweep  oscillator  is  then  set  at  the  assigned  frequency  by  the 
computer  and  the  power  output  level  is  set  to  its  lowest  possible  point.  The 
waveguide  wavelength  is  calculated  for  the  assigned  frequency  and  a  search 
process  is  begun  to  find  the  standing  wave  minimum  nearest  the  sample.  The 
search  is  a  four-part  iteration  process. 

With  the  probe  positioned  at  its  initial  location,  the  microwave  power  is 
set  very  low,  to  less  than  70  dBm.  Each  iteration  commands  the  source  to 
increase  the  power  of  the  microwave  signal  in  increments  while  the  SWR  meter 
outputs  a  voltage  to  the  multimeter.  When  the  multimeter  detects  a  reading  of 
100  mV  or  more,  the  power  level  of  the  source  is  recorded  and  then  reset  low, 
again  less  than  70  dBm.  The  stepper  motor  is  then  commanded  to  move,  the  probe 
one  increment  (a  step).  The  powerjs  again  increased  until  the  100  mV  can  be 
read  from  the  multimeter,  the  power  level  is  recorded  and  the  power  is  set 
low,  and  the  probe  is  again  moved.  The  third  power  reading  is  compared  to  the 
first  two  to  see  how  the  power  level  has  changed  before  the  probe  is  moved 
again.  The  program  looks  for  the  position  where  the  power  level  decreases 
after  two  consecutive  rises  in  recorded  power  level.  The  highest  power  level 
reading  is  assumed  to  be  the  estimate  of  SWR  minimum  for  that  iteration. 

In  the  first  iteration,  the  wavelength  is  divided  Into  10°  step  sizes. 
Since  there  are  two  standing  waves  for  each  waveguide  wavelength  or  360°,  at 
most  twenty  steps  would  be  necessary  for  the  first  iteration.  This  would  only 
be  necessary  in  the  case  that  the  probe  started  at  or  less  than  20°  to  the 
left  of  a  standing  wave  minimum,  thereby  requiring  a  search  over  200°.  The 
probe  is  placed  half  of  the  initial  step  (5°)  before  the  maximum  power  reading 
location  -  the  first  estimate  of  the  standing  wave  minimum  from  the  first 
iteration.  A  new  step  size  is  calculated  for  steps  of  1°,  and  again,  after 
the  power  has  increased  for  two  steps  and  then  drops  off,  the  probe  is  located 
a  half  step  (.5°)  before  the  maximum  power  reading.  A  new  step  size  of  .1°  is 
calculated  and  used  in  the  next  iteration  using  the  same  process  as  above.  A 
final  step  size  of  .05°  is  used  in  the  last  iteration  to  locate  the  best 
estimate  of  the  minimum  position.  The  probe  is  then  located  at  that  position. 

A  new  process  is  begun  to  locate  the  "3  dBm  locations."  These  are  the 
points  on  either  side  of  the  standing  wave  minimum  where  the  SWR  meter  power 
reading  would  be  3  dBm  greater  than  at  the  SWR  minimum.  The  program  finds 
these  points  by  turning  the  power  level  of  the  source  down  3  dBm  from  the 
maximum  power  reading  at  the  standing  wave  minimum.  The  probe  is  moved  away 
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from  the  minimum,  both  positively  and  negatively,  using  step  sizes  of  .025° 
until  the  3  dBm  positions  are  located.  This  occurs  on  either  side  of  the 
minimum  position  where  the  voltage  reading  returns  to  the  level  obtained  when 
the  probe  was  located  at  the  standing  wave  minimum.  This  completes  the 
baseline  measurement  without  the  sample  in  place.  The  probe  carriage  is  then 
returned  to  its  original  starting  location. 

The  computer  at  this  point  prompts  the  user  to  load  the  sample  in  place 
in  the  sample  holder.  If  a  porous  sample  is  being  tested,  a  known  weight  of 
material  is  loaded  and  then  pressed,  if  required,  to  achieve  a  higher 
density.  After  the  sample  in  its  sample  holder  is  attached  to  the  circuit, 
the  process  is  begun  anew  to  locate  the  first  new  standing  wave  minimum  to  the 
left  of  the  starting  position.  The  locations  of  the  "3  dBm  points",  both  with 
and  without  the  sample  in  place,  are  sent  to  the  portion  of  the  program  that 
calculates  the  dielectric  properties  (lines  1870  to  3170  of  program  shown  in 
Appendix  A).  The  original  program,  as  adapted  by  R.  Watkins  of  this  Center, 
is  shown  in  Appendix  B.  That  program  takes  manually  input  values,  previously 
determined  from  the  adjustment  of  the  probe  carriage,  to  calculate  the 
dielectric  parameters.  See  Reference  2  if  further  details  are  needed 
regarding  this  program. 

When  the  calculation  of  the  dielectric  properties  is  complete,  the  first 
three  calculations  of  possible  property  values  are  printed  out.  The  computer 
then  asks  if  more  values  are  desired;  if  so  press  "Y",  if  not  press  "N"  and 
the  program  will  stop.  If  an  estimate  of  the  dielectric  constant  is  not 
known,  then  the  process  will  need  to  be  repeated  with  a  different  sample 
length.  The  intersection  of  the  two  data  sets  is  the  correct  dielectric 
constant  value. 


NSWC  TR  86-46 


CHAPTER  4 
FUTURE  WORK 

The  frequency  source  plug-in  used  for  early  studies  had  a  frequency  range 
of  0.1  to  8.4  GHz  and,  therefore,  only  had  a  usable  range  of  8.2  to  8.4  GHz  in 
X-band  waveguide.  A  C-band  setup  could  be  developed  to  better  utilize  the 
frequency  capabilities  of  the  available  plug-in.  For  the  C-band  system,  the 
program  would  be  adapted  to  make  automated  measurements  over  the  entire  C-band 
range  of  frequencies  at  specific  intervals.  Then  the  variation  of  dielectric 
properties  over  a  range  of  frequencies  for  a  specific  sample  could  readily  be 
seen.  Recently,  a  frequency  source  plug-in  became  available,  covering  the 
entire  X-band  range.  For  this  plug-in,  the  program  will  be  changed  to  make 
measurements  automatically  over  the  entire  X-band  range  at  specific  frequency 
intervals.  Then,  with  porous  samples,  dielectric  constants  can  be  determined 
at  various  frequencies  and  sample  densities. 

Temperature  could  also  be  introduced  as  an  additional  variable.  The 
sample  holder  could  be  placed  into  an  oven  arrangement  to  allow  the 
temperature  to  be  increased  incrementally  as  required  by  the  program.  The 
temperature  measurement  of  the  sample  could  be  accomplished  with  the  use  of  a 
thermocouple  to  measure  the  temperature  and  an  IEEE-488  controllable  power 
supply  to  control  the  temperature  of  the  oven.  Both  of  these  could  be  easily 
added  to  the  existing  program.  It  would  be  possible  at  that  time  to  acquire 
data  as  the  frequency  varies  in  conjunction  with  changes  in  temperature  for 
each  specific  sample  density  under  test. 

The  original  developmental  use  for  this  setup  was  to  tailor  a  material 
with  equal  permittivity  and  permeability  for  dielectric  antennas;  this 
research  will  continue.  It  is  also  possible  for  the  magnetic  properties  of 
the  materials  to  be  calculated  using  the  type  of  procedure  used  here.  This 
requires  a  different  equation  to  be  solved  in  the  calculation  of  magnetic 
permeability  and  dielectric  constant,  and  the  measurement  of  the  sample  at  two 
different  positions  in  the  waveguide. 


CHAPTER  5 


VARIATIONS  OF  DIELECTRIC  PROPERTIES  WITH 
VARIATIONS  OF  TEMPERATURE  AND  HUMIDITY 


A  porous  sample  is  considered  to  be  a  mixture  of  solid  material  and  air 
where  the  solid  material  is  diluted  from  its  theoretical  maximum  density 
(TMD).  Its  electrical  permeability  or  relating  dielectric  constant  is  then 
calculated  using  the  Lichtenecker-Rother  and  the  Landau-Lifshitz-Looyenga 
formulas  for  combining  materials  of  variable  volume  percentages. 

No  adequate  theory  is  available  to  predict  the  effect  of  temperature 
change  on  the  dielectric  constant,  and  measured  materials  have  exhibited  both 
positive  and  negative  temperature  coefficients.  Fortunately,  at  microwave 
frequencies,  most  of  the  molecular  resonances  have  occurred  and  the 
temperature  effects  are  small.  Changes  of  less  than  .01  percent/°C  in  the 
dielectric  constant  are  normally  expected,  however,  temperature  coefficients 
should  be  measured  for  each  new  material  used. 

The  effects  of  humidity  or  absorbed  water  in  the  sample  represent  a 
potentially  large  variation  in  dielectric  constant  because  water  has  a  high 
relative  dielectric  constant  ("55  at  10  GHz).  To  add  to  this,  its  temperature 
coefficient4  near  room  temperature  is  about  1  percent/°C  and  its  frequency 
variation4  is  about  3  percent/GHz.  Fortunately,  in  a  sealed  waveguide  system, 
such  as  is  used  in  the  experiments,  the  external  environment  is  of  little 
concern.  Our  samples  were  prepared  in  a  room  with  a  controlled  temperature 
and  low  relative  humidity.  For  these  reasons,  the  effects  of  temperature  and 
humidity  have  been  neglected  in  our  analysis,  though  are  reported  for  future 
considerations. 


CHAPTER  6 


RESULTS  AND  DISCUSSION 


Although  this  technique  gives  good  values  for  the  real  part  of  the 
dielectric  constant,  it  is  known  to  be  less  accurate  for  the  determination  of 
values  for  loss  factor  or  loss  tangent.  To  obtain  more  accurate  values  for 
the  loss  parameters,  other  methods  of  testing  could  be  used,  but  they,  in  turn 
are  less  accurate  in  determining  the  real  values  of  dielectric  constant.  Some 
data  printouts  are  shown  in  Appendix  C  as  an  illustration  of  the  performance 
of  this  program.  The  first  two  examples  were  run  to  show  that  this  procedure 
would  return  a  dielectric  constant  value  of  1.0  for  an  air  filled  sample 
holder.  For  the  two  sample  holder  lengths  used,  the  dielectric  constant  data 
intersects  at  the  value  of  1.000  and  the  accepted  value  for  air  is  in  fact 
1.00064. 

The  next  example  of  data  printout  is  for  a  sample  of  ceramic  foam  with  a 
dielectric  constant  known  to  be  near  1.5.  Our  procedure  confirmed  this  by 
-  calculating  the  most  likely  value  of  dielectric  constant  to  be  1.507.  A 
further  demonstration  of  the  effectiveness  of  this  procedure  to  find  the 
correct  value  of  dielectric  constant  was  shown  when  a  sample  of  polysulfone 
was  tested.  Polysulfone  has  been  extensively  measured  by  the  Massachusetts 
Institute  of  Technology  (MIT),  Laboratory  for  Insulation  Research  and  their 
findings  show  the  variation  of  dielectric  constant  with  temperature  and 
incident  microwave  frequency.  The  temperature  and  frequency  parameters  which 
most  closely  matched  our  test  parameters  occurred  at  25°C  (77°F)  and  8.515  GHz. 
The  MIT  value  for  dielectric  constant  was  given  as  2.980.  Our  procedure 
determined  the  dielectric  constant  at  a  temperature  of  79°F  and  an  incident 
microwave  frequency  of  8.40  GHz  to  be  2.971. 

The  automated  testing  apparatus  was  used  as  a  research  tool  to  determine 
the  dielectric  properties  of  a  porous  inert  material.  The  material  chosen  for 
investigation  was  Teflon  7C.  Samples  were  tested  at  a  constant  frequency  of 
8.30  GHz  in  a  well  controlled  environment,  where  temperature  and  humidity  were 
kept  constant.  A  variety  of  loading  densities  were  tried  to  give  the 
variation  of  the  dielectric  constant  with  loading  density.  This  data  for 
Teflon  7C  is  plotted  in  Figure  2.  It  was  observed  from  the  plot  that  the  data 
follows  the  empirical  exponential  formula 

v  v 

sr  *  n  e2  2 

where  er  *  relative  dielectric  constant 

e,  *  dielectric  constant  of  the  material  at  TMO 
vj  *  solid  fraction  of  the  material.  Percent  TMD/100 
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FIGURE  2.  VARIATION  OF  DIELECTRIC  CONSTANT  WITH  PERCENT  TMD  FOR 
TEFLON  7C  POWDER  AT  8.3  GHz 


e2  »  dielectric  constant  of  mixing  material,  air,  in  this  simple  case 
v2  “  *  “  V1  *  volume  fraction  of  mixing  material 

This  formula  was  derived  previously  by  Lichtenecker-Rother^  in  the  early 
1930 's  to  describe  the  dielectric  constant  resulting  from  the  mixture  of 
various  components.  In  the  simple  case  where  e2  an<1  v2  are  for  air,  Its 
dielectric  constant  is  1.0,  therefore  e2  raised  to  the  v2  power  is  still  one, 
and  can  be  neglected. 

The  apparatus  again  was  used  to  evaluate  the  same  dielectric  constant 
variation  for  Class  D  HMX.  The  data  for  this  testing  Is  shown  in  Figure  3. 
Since  most  of  our  testing  is  done  with  Class  0  HMX  porous  beds  initially 
loaded  at  73  percent  theoretical  maximum  density  (percent  TMD),  a  large  number 
of  trials  were  run  at  that  density  to  obtain  an  average  dielectric  constant 
for  the  frequency  of  8.30  GHz.  After  determining  that  the  empirical 
exponential  equation  fit  the  data  well,  the  number  of  tests  run  on  any 
particular  density  need  only  be  one,  assuming  a  range  of  densities  are  tested. 

The  given  formula  and  the  automated  dielectric  constant  measurement 
apparatus  were  again  used  to  describe  the  variation  of  dielectric  constant 
with  density  for  various  mixtures  of  melamine  (inert)  and  aluminized 
powders.  The  mixtures  contained  50u  particle  size  melamine.  Measurements 
were  made  of  melamine  alone,  for  use  as  a  base  line,  and  for  10,  20  and  30 
percent  of  aluminum  by  weight  in  the  melamine.  The  dielectric  constant  of 
aluminum,  a  conductor,  is  not  readily  known,  but  the  dielectric  constant  of 
the  mixture  can  be  easily  measured,  thus  allowing  the  use  of  the  lichtenecker- 
Rother  formula,  only  now  the  variables  are  slightly  changed.  Now,  e,  is  the 
dielectric  constant  at  TMD  for  the  mixture  of  melamine  and  aluminum  as 
determined  empirically  from  experimental  data,  e_  is  the  relative  dielectric 
constant  of  the  mixture  and  v^  is  now  the  solid  fraction  of  the  entire 
mixture.  This  data  Is  shown  graphically  in  Figure  4.  All  data  for  Figures  2, 
3,  and  4  are  given  in  Table  1. 

Once  the  value  of  relative  dielectric  constant  at  TMD  for  each  mixture 
has  been  computed,  these  TMD  values  can  be  replotted  as  the  dielectric 
constant  versus  the  percentage  of  aluminum  from  0  to  100  percent.  The  values 
at  low  percentages  can  be  reinserted  into  the  mixture  formula  to  calculate  the 
dielectric  constant  of  100  percent  aluminum  or  a  curve  can  be  fitted  through 
the  points  and  extrapolated  to  100  percent.  Either  method  yields  a  value  of 
the  relative  dielectric  constant  of  aluminum  to  be  approximately  12.0.  This 
is  the  first  reporting  of  the  experimental  value  of  the  dielectric  constant  of 
such  a  highly  conductive  material  known  to  the  author.  Normally,  the 
conductivity  and  associated  losses  of  aluminum  or  other  conductors  preclude 
such  a  measurement.  Not  only  can  this  mixture  method  be  used  for  the 
measurement  of  metals  and  semiconductors,  it  can  also  be  used  to  measure  any 
other  material  with  a  very  high  dielectric  constant. 

This  automatic  apparatus  for  the  measurement  of  the  dielectric  constant 
has  given  us  a  necessary  tool,  enabling  the  measurement  of  a  wide  variety  of 
porous  and  solid  materials,  to  be  practically  done  in  a  modest  amount  of 
time.  Further,  this  has  allowed  the  use  of  impedance  matching  techniques 
which  had  previously  been  used  only  in  broad  band  ranges. 
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FIGURE  3.  VARIATION  OF  DIELECTRIC  CONSTANT  WITH  PERCENT  TMD  FOR 
CLASS  0  HMX  AT  8.3  GHz 
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FIGURE  4.  VARIATION  OF  DIELECTRIC  CONSTANT  WITH  PERCENT  TMD  FOR 
VARIOUS  MIXTURES  OF  MELAMINE  AND  ALUMINUM  AT  10GHz 


TABLE  1.  EXPERIMENTAL  DIELECTRIC  CONSTANT  DATA 
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TABLE  1.  (Cont.)  EXPERIMENTAL  DIELECTRIC  CONSTANT  DATA  (Cont) 

Dielectric  Relative 

Sa*P'e  Length  (n»)  Density  (g/cc)  %  TMD  Constant  Temp  (F)  Humidity  ( 

Melamine  12.14  0.791  46.1  2.307  75  62 
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APPENDIX  A 


AUTOMATED  DIELECTRIC  CONSTANT  MEASUREMENT  PROGRAM 
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THIS  PROGRAM  IS  STORED  AS  "AUTO-E" 

ON  THE  TAPE  LABELED  “BC  85“ 

APRIL  1985,  Brian  Clancy  and  A1  Krai) 

THE  SAMPLE  IS  INSERTED  IN  THE  END  OF  A  SECTION  OF  WAVEGUIDE  AFTER 
THE  PROGRAM  IS  BEGUN,  WHEN  IT  IS  INDICATED  THAT  YOU  DO  SO. 

A  SLOT  IN  THE  TRANSMISSION  LINE  ALLOWS  THE  VSUR  AND  MINIMUM 
POSITION  OF  THE  RESULTING  STANDING  WAVE  TO  BE  MEASURED. 

THIS  PROGRAM  IS  FOR  AN  HP  9845B  COMPUTER  AND  THE  INSTRUMENTATION  FOR 
AN  EXPERIMENT  TO  DETERMINE  THE  DIELECTRIC  PARAMETERS  OF  ANY  MATERIAL. 
THIS  PROGRAM  IS  DESIGNED  TO  KEEP  THE  POWER  LOW  SO  THAT  THE  SLOTTED 
LINE'S  PROBE  CRYSTAL  IS  ALWAYS  OPERATING  IN  ITS  SQUARE  LAW  REGION. 

THE  PROGRAM  ADJUSTS  THE  POWER  TO  KEEP  THE  VOLTAGE  LOW,  THEREFORE, 

THE  PROBE  WILL  BE  LOCATED  AT  A  MINIMUM  WHEN  THE  POWER  IS  HIGHEST. 

ADDITIONAL  INSTRUMENTATION: 

Klinger  Stepper  Motor  Controller,  CC-1  and  Mi cro-Control *  Motor 
HP  83S0A  Sweep  Oscillator  with  83S25A  RF  plug-in,  .01-8.4  GHz 
HP  .3438A  Digital  Multimeter  connected  to  an  HP  415D  SWR  Meter 
Appropriate  Slotted  line  probe  carriage  and  Associated 
Waveguide  for  the  Frequency  Range  of  Interest 


IF  THE  TRANSMISSION  LINE  IS  CHANGED  FROM  RECTANGULAR  X-BAND,  CHECK 
TO  MAKE  SURE  THE  PROGRAM  LINES  CONTAINING  RECTANGULAR  X-BAND  DATA 
HAVE  BEEN  CHANGED  TO  REFLECT  THE  NEW  TRANSMISSION  LINE  SI2E  and/or 
WAVE  VELOCITY! 

IF  THE  PROGRAM  IS  NOT  CONVERGING  TO  A  POWER  MAXIMUM,  PRESS  STOP  ON 
.COMPUTER,  PRESS  COUNTER  RESET  ON  THE'  KLINGER,  THEN  RUN  ON  COMPUTER. 
THIS  WILL  RESTART  THE  PROGRAM  AT  A. DIFFERENT  POSITION 


PRINTER  IS  16  !  PRINTER  16  IS  THE  CRT  SCREEN  OF  THE  9845B 

OPTION  BASE  1 

DIM  Dba<40) , Vo! ts<40) ,Pos<40> , Sampl SI80I,Xx(9,2),Dc(9), £pp(9) , TadC9> 

DIM  F<3, 30) 

DEF  FNSi nh(X)«. 5*<EXP<X)-£XP<-X) )  !  USED  IN  DIELECTRIC  CONST.  CALC. 
DEF  FNCosh(X)«.5*<EXP<X)*EXP(-X))  !  USED  IN  DIELECTRIC  CONST.  CALC. 
Al-22.86  !  Width  of  Rectangular  Waveguide  in  ••  (X  Band) 

81*10. 16  !  Height  of  Rectangular  Waveguide  in  aa  (X  Band) 

Slot-3. 2  !  Rectangular  Slotted  line  slot  width  (•■)  (X  Band 

OUTPUT  9; “Request  time"  !  FROM  THE  98035A  REAL  TIME  CLOCK 
ENTER  9; Month, Day, Hour, Mi nute, Second 

GOSUB  3910  !  SENDS  US  TO  A  SUBROUTINE  TO  HELP  US  PRINT  OUT  DATE  1  TIME 

BEEP  ‘BELL 

PRINT  "IS  THE  SLOTTED  SECTION'S  PROBE  NEAR  THE  RIGHT  SIDE  OF  CARRIAGE?" 
PRINT  “  RESET  THE  KLINGER  BY  PRESSING  THE  COUNTER  RESET  BUTTON- 

INPUT  “What  Year  is  it  now?  Enter  only  the  last  two  digits", Year 
INPUT  "Room  Temperature  is  what  at  this  time?  (Deg.  F)",Teap 
INPUT  "Relative  Huaidity  in  the  room  is  what  at  this  time?  CS)",Rh 
PRINT  CHRSC12)  !  FORM  FEED  TO  CLEAR  THE  CRT  SCREEN 
Count*l 

IF  Count-1  THEN  GOTO  740 

INPUT  "Is  the  sample  in  place  on  the  set-up?  (Y  or  N>  CAPS*,Sa»p* 

IF  Samp*-" Y"  THEN  GOTO  630 
BEEP  !  BELL 

PRINTER  IS  16 

PRINT  “PUT  SAMPLE  IN  PLACE.  NOW !  " 

PRINTER  IS  0 
GOTO  550 

BEEP  !  BELL 

PRINTER  IS  16 

PRINT  CHRSC12)  !  FORM  FEED  TO  CLEAR  THE  CRT  SCREEN 
PRINTER  IS  0 
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1858 
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1898 
1188 
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1128 
1138 
1148 
1158 
1158 
1178 
1188 
1198 
1288 
1218 
1228 
1238 
1248 
1258 
1268 
1278 
1288 
1298 
1388 
1318 
1328 


SAMPLE  DESIGNATION  <  1  ess  than  38  Ch4r*et*ri)a,SMptt 


INPUT 
PRINT 
PRINT  Sup  1  * 

INPUT  "LENGTH  OF  SAMPLE  <mm)“ ,Li 
PRINT  “Length  of  Sample  is: 

PRINT 
COTO  768 

INPUT  “FREQUENCY  <GHs)",F 
OUTPUT  7, 19; *CW“,F, "GZ" 

Dbm— 40. 0 

OUTPUT  7, 19; "PL"f  Dbm, "DM" 

EOL  7; CHRSC 13)  !  ADDRESS 

!  THE  EOL  FUNCTION  DEFINES 


;  Le;  “<»(•>• 


»  DEFINES  THE  FREQUENCY  TO  BE  USED 


!  SETS  THE  FREQUENCY  ON  THE  HP  S350A  SOURCE 
!  DEFINE  THE  dBm  LEVEL  TO  BE  -48.8  dBm 
!  SETS  dBm  LEVEL  ON  THE  3358A  SOURCE 
7  IS  THE  KLINGER  STEPPER  MOTOR  CONTROLLER 
A  NEW  END  OF  LINE  CHARACTER,  IN  THIS  CASE 
!  CHRS (13)  DEFINES  A  CARRIAGE  RETURN  WITH  NO  LINE  FEED,  AS  NECESSARY 
!  FOR  THE  KLINGER  STEPPER  MOTOR  CONTROLLER. 

!  THE  NEW  END  OF  LINE  STRING  IS  STORED  AS  IMAGE  "L*. 

!  THE  "USING  "*, K, L" ,  IS  USED  TO  DEFINE  THE  USE  OF  THE  NEW  EOL  STRING, 

!  THE  "*"  DISABLES  THE  CARRIAGE  RETURN  AND  LINE  FEED  COLD  EOL  STRING) , 

!  THE  "L“  ALLOWS  US  TO  USE  OUR  EOL  STRING  STORED  AS  IMAGE  "L". 

Laad40-299.6966/F  !  2.997E8  M/'SEC  IS  THE  SPEED  OF  E-M  RAD.  IN  AIR  9  ST P 

!  THIS  LAST  EQUATION  CALCULATES  THE  FREE  AIR  WAVELENGTH  OF  ANY 
!  FREQUENCY  OF  ELECTROMAGNETIC  RADIATION. 

Lamdac-2*A1  ■  THE  CUTOFF  FREQ.  FOR  RECTANGULAR  WAVEGUIDE  CTE18  MODE) 

Lud4g<Lud»e/(l*(Lud*8/LudtC)A2)".S  !  WAVEGUIDE  WAVELENGTH 

PRINTER  IS  8  !  PRINTER  8  IS  THE  HARD  COPY  PRINTER  OF  THE  9845B 

IF  Count-2  THEN  GOTO  1878 

PRINT 

PRINT 

IMAGE  "Date:  ", 22, A, 22, A, 22 

PRINT  USING  958; Month, SI aahS, Day, SI ash$, Year 
IMAGE  "Tim*:  ",22,8,22,8,22,"  At  atart  of  teat" 

PRINT  USING  978; Hour , Col *, Mi  nut*, Cot *, Second 
PRINT 

PRINT  "Rel at i ue  Huni di ty-" ; Rh; "S* , "Temperature-" ; Temp; "Dog.  F* 

PRINT 

IMAGE  "Frequency  selected  for  this  test  i s" , DDD. ODD, "  CCHz)"' 

PRINT  USING  1828; F 

IMAGE  "Waveguide  Wavelength  -" , DDD. DDD, " " 

PRINT  USING  1048;  l^amdag 
PRINT 

Step-INTCLaadag^36«1888)  !  READINGS  TAKEN  EVERY  18  DEGREES  OF  WAVELENGTH 


GOSUB  3788 
Latep-Step 
Pm*-**" 

Datepl-2 
Dstep2«l 
Ds t  ep3- . 7 
Dbl ess-28. 8 
GOSUB  4848 
Step-INT<Lstep/2) 

GOSUB  3788 
Dbm-Dbm-7 

OUTPUT  7, 19; “PL  -28  DM" 
GOSUB  3688 
Pm*-"+" 

Dstepl-l 
Dstep2«. 5 
Dstep3>. 1 
Dbless-18 

St ep-INT(La»dag/360*l080> 

Lstep-Step 

Dbm-Dbm- 1 8 

OUTPUT  7, 19; "PL  -5  DM" 
GOSUB  3708 
GOSUB  4840 
Step«lNT<L*tep'2) 


!  DEFINE  THE  MOVE  COMMAND  FOR  THE  KLINGER  ‘N  Step" 
Lstep  -  LAST  STEP  SUE 
SET  DIRECTION  OF  NEXT  MOVE 

THESE  THREE  Dstep  VALUES  ARE  USED  TO  ADJUST  THE  8350A 
dBm  LEVEL  IN  THE  SUBROUTINE  THAT  ZEROS  IN  ON  dBm 


SUBROUTINE  THAT  FINDS  THE  NEXT  ESTIMATE  OF  MINIMUM 
DEFINE  THE  MOVE  COMMAND  FOR  THE  KLINGER  "N  Step" 


MOVE  PROBE  TO  POSITION  1/2  STEP  BEFORE  MINIMUM  ESTIMATE 


EACH  STEP  IS  1  DEGREE 


DEFINE  THE  MOVE  COMMAND  FOR  THE  KLINGER  "N  Step" 
FIND  NEXT  ESTIMATE  OF  MINIMUM  POSITION 


A-3 


I»i 

\ 


1330 

1340 

13S0 

1360 

1370 

1300 

1390 

1400 

1410 

1420 

1430 

1440 

1430 

1460 

1470 

1480 

1490 

1300 

1310 

1320 

1330 

1340 

1330 

1360 

1570 

1380 

1590 

1600 

1610 

1620 

1630 

1640 

1630 

1660 

1670 

1680 

1690 

1700 

1710 

1720 

1730 

1748 

1730 

1760 

1770 

1780 

1790 

1300 

1810 

1320 

1830 

1840 

1330 

1860 

1879 

1338 

1890 

1900 

1919 

1920 
1930 
1940 
1930 
I960 
1970 
1990 
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COSUB  3680  !  DEFINE  THE  MOVE  COMMAND  FOR  THE  KLINGER  "N  Step” 

GOSUB  3600  !  MOVE  THE  PROBE  TO  A  POSITION  i/2  STEP  BEFORE  MIN.  EST. 

Pm»m"f“  !  SET  DIRECTION  OF  NEXT  MOVE 

Dstepl*. 5 
Dst  ep2* •  1 
Dst ep3*. 92 
Dbl ess*l 

Step*INT<L**d*g/3600*1000>  !  EACH  STEP  *  .1  DEGREE 

GOSUB  3790  !  DEFINES  THE  MOVE  COMMAND  FOR  THE  KLINGER  "N  Stef 

Lstep»Step 
Dbm*Db«-l 0 

OUTPUT  7, 19; "PL  3.0  DM" 

GOSUB  4040  !  FIND  THE  NEXT  ESTIMATE  OF  MINIMUM  POSITION 

Step*Step*2 

GOSUB  3700  !  DEFINES  THE  MOVE  COMMAND  FOR  THE  KLINGER  *N  St*| 

GOSUB  3600  !  MOVE  TO  POSITION  2  STEPS  BEFORE  MINIMUM  ESTIMATE 

t  set  DIRECTION  OF  NEXT  MOVE 

Dstepl*. 2 
Dstep2*. 08 
Dstep3*. 04 
Dbl ess*. 5 

St ep*I NT (Lam dag/ 7200* 1008 )  !  EACH  STEP  IS  .05  DEGREE 

IF  Step*0  THEN  Step*l 

GOSUB  3780  !  DEFINES  THE  MOVE  COMMAND  FOR  THE  KLINGER  “N  St*| 

Lstep*Step 
Dbm*Dbm-l 

OUTPUT  7,  19; *PL",Dbm, "DM* 

GOSUB  4040  !  FIND  THE  NEXT  ESTIMATE  OF  MINIMUM 

Dbm-Dbm-3  !  SET  THE  dBm  POWER  DOUN  3  dBm  FROM  L.AST  READING 

OUTPUT  7, 19; “PL", Dbm, "DM" 

WHIT  1000  !  WAIT  l  SECOND  UNTIL  VOLTMETER  HAS  Tii"*?  TO  SETTLE 

!  Lst ep*St *p 
Step*INT<Lstep/2> 

Lstep*Step 

IF  Step*0  THEN  Step-l 


•N  Step" 


!  DEFINES  THE  MOVE  COMMAND  FOR  THE  KLINGER  *N  Step* 
!  MOVE  TO  POSITION  2  STEPS  BEFORE  MINIMUM  ESTIMATE 
!  SET  DIRECTION  OF  NEXT  MOVE 


!  EACH  STEP  IS  .03  DEGREE 


*N  Step* 


GOSUB  3700 
GOSUB  3730 
GOSUB  3860 
Pt3dbl*Pos 
Step*I 
COSUB  3700 
P»*«"f" 
GOSUB  3600 
Step*Lstep 
WAIT  2090 
GOSUB  3780 
GOSUB  3730 
GOSUB  3360 
Pt3db2*Pos 


!  DEFINES  THE  MOVE  COMMAND  FOR  THE  KLINGER  "N  S 
!  MOVE  KLINGER  TO  FIND  THE  3  dBm  POSITION 
!  ASK  THE  KLINGER  THE  POSITION  OF  THE  PROBE 
!  POINT  OF  FIRST  3  dBm  POSITION 

!  FIGURE  STEP  TO  MOVE  BACK  TO  THIRD  MINIMUM  EST. 
!  DEFINE  THE  MOVE  COMMAND  FOR  THE  KLINGER 
!  SET  DIRECTION  OF  NEXT  MOVE 
!  MOVE  PROBE  BACK  TO  THIRD  EST.  OF  MINIMUM 


N  Step" 


!  WAIT  2  SECONDS  BEFORE  SEARCHING  FOR  2nd  3  dBm  POSITION 
!  DEFINE  THE  MOVE  COMMAND  FOR  THE  KLINGER  “N  Step" 

!  MOVE  KLINGER  TO  FIND  THE  3  dBm  POSITION 
•  ASK  THE  KLINGER  THE  POSITION  OF  THE  PROBE 
!  POINT  OF  SECOND  3  dBm  POSITION 


OUTPUT  7, 19; "PL  -72  DM" 
St  ep*Pos 
GOSUB  3700 
PmS*“ - “ 

GOSUB  3690 


!  SET  POWER  LEVEL  OF  3350A  BACK  LOW 


!  THIS  PORTION  OF  THE  PROGRAM  IS  BASED  ON  THE  NOV.  1972  REPORT 
!  ARS-NC-4  FROM  THE  AGRICULTURAL  RESEARCH  SERVICE.  IT  WILL  FIND 
!  DIELECTRIC  CONSTANT,  LOSS  TANGENT,  AND  LOSS  FACTOR. 


IF  Count  *2  THEN  GOTO  1980 


LZ**Pt3dbl/1000  !  Lz*  IS  ONE  3  dBm  POSITION  U/O  SAMPLE 

Rz**Pt3db2'l000  !  Rz*  IS  ONE  3  dBm  POSITION  U/O  SAMPLE 

PRINT  "Position*  ot'  3  dBm  pts  w/o  sample  ere  " ;  Lz*;  Rz*.;  "mm" 
Count  *Count  *  1 


GOTO  330 
Rzs-Pt3dbl '10O0 


!  Rzs  IS  ONE  3  dBm  POSITION  WITH  SAMPLE  IN  PLACE 


ft 


fti 


vv.'-.v-.v 


m 


1990  Lzs*Pt  3db2'"  1000  !  Lzs  13  ONE  3  dBm  POSITION  WITH  SAMPLE  IN  PLACE 

2000  PRINT  "Position*  of  3  dBm  pts  with  s*mp1«  *r* " ; Rzs ; Lz* ; “ mm “ 

2010  Sd“3 
2020  Dt>*3 

2030  Wl—  2 

2040  W 1 «U 1 ♦2  '  ASSIGN  AN  INITIAL  VALUE  FOR  SHIFT 

20S0  n«x*0 

2060  Sc«l*<Slot*L*md*g>''2/<8*PI*Bl*Al''3>  !  SLOTTED  LINE  WAVELENGTH  CORRECTION 

2070  2*-.3*<Lz*+Rz*>  !  MINIMUM  POSITION  WITHOUT  SAMPLE 

2080  2l»Z*  !  REDEFINE  MINIMUM 

2090  Q-l 

2108  IF  A8S(Q*(L*«d*g'2>-<2*+LAmd*g/'2)  ><L*»d*g/4  THEN  2130 
2110  Q*Q*  1 
2120  GOTO  2100 

2130  2*«Q*<L*md*g^2>  !  MIN  POS  W/O  SAMPLE  IS  INTEGRAL  *  L*md*g/2 ' S 

2140  2*«.5*<Lz*->-Rz*>  !  MINIMUM  POSITION  WITH  SAMPLE  IN  PLACE 

2130  Dz*»ABSCLz*-Rz»>>Sc  !  CORRECTED  SEPERATION  BTW  3dB  PTS  U'O  SAMPLE 

2160  Dz*-ABS(Lzs-Rzs)''Sc  !  CORRECTED  SEPERATION  BTW  3dB  PTS  WITH  SAMPLE 

2170  St«<21-2*>'Sc  !  CORRECTED  SEPERATION  BTW  MINIMUMS 

2180  Swr«< 1. 9932623-COS <P I *Dz*/L*md*g> *2 3/S  INC  PI *Dz*/'L*«d*g> 

2190  Arg-(.9932623/'(Sur''2-l>>''.5 

2200  Dz*mL*«nd*g*ASN<  Arg>^PI 
2210  R»<2*-L*-St)/2* 

2220  Ppr*PI*<Dz*-R*Dz*)/L**d*g 

2230  S«r-( 1 . 99S2623-C0S<Ppp>~2>~. 5/SlN(Ppr > 

2240  Ph  i  *  1  /Sur* 

2230  N-0 

2260  IF  -<L«-*-St)  +  .5*L*nd»g*N>0  THEN  2290 

2270  N-N+i 

2280  GOTO  2260 

2290  20— <St*L*)  +  .  5*L*»d*g*N 

2308  U**2*PI*Z0/'L*»d*g 

2310  A-TAN<U>+£Phi-2-l> 

2320  B«Phi*<  1+TPN<U>'"2J 

2330  C*<A~2+B‘'2>-.5*<L*«d*g^<2*PI*L*-iKl'f(Phi  ♦TAN<U>>-'2>>  > 

2340  C»C*<A/ABS(A)) 

2330  IF  C<1  THEN  K-2 

2360  FOR  Y-K  +  1+W1*2  TO  6+K+Wl*2  STEP  2 

2378  Ml-Y*PI^2 

2380  M2-(Y-2>*PI^2 

2390  IF  M2< 0  THEN  M2-0 

2400  X2«(Ml-M2)/'2+M2 

2410  FU«TAN(X2VX2-C 

2420  IF  Fu>0  THEN  M1-X2 

2430  IF  Fu< “0  THEN  M2-X2 

2440  IF  Ml-M2(. 000000001  THEN  2460 

2430  GOTO  2400 

2460  M«M+t 

2470  Xx(M,l)«X2 

2480  Ad»Xx(M,  t)*L*md*g*Phi*<l*TAN(U)''2) 

2490  D*n«2*P I 4L**< l+TAN(Xx(M, 1 > > -2- TAN ( Xx ( M , 1 ) >^Xx(M, 1 > ) 

2300  XxCM, 2>*Ad'D*n 

2310  NEXT  Y 

2320  Q*L*md».g^  ( 2*P  I  *L* ) 

2330  D*n«Q/< l*<Phi *TAN(U> >~2> 

2340  U0*A*D*n 

2330  V0— B*D*n 

2360  FOR  X»l  TO  3 

2370  J-0 

2380  X6»Xx  <  X , 2  > 

2398  X7*Xx(X,  1  'j 

2600  D-.01 

2610  GOSUB  3930 

2620  Di  t"(  (U0-U*;''2*'.V0-V*;''2)^.  3 
2630  IF  D i 1 < . 00001  THEN  2920 
2640  K«0 
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2630  X6*Xx  <  X , 2 / *D 

2660  X7*Xx  <  X ,  1  ) 

2570  GOSUB  3930 

2680  Di2»<<U0-U»)^2+<VO-V*>-2>''.  5 

2690  K«K+l 

2780  IF  Di2>«Di 1  THEN  2740 

2710  Xx<X,2)*X6 

2720  Di l*Di  2 

2730  GOTO  2630 

2740  IF  K>  1  THEN  2770 

2730  X6«Xx<X,2>-D 

2760  GOTO  2670 

2770  K«0 

2780  X7-Xx<X, 1)+D 

2790  X6«Xx<X,2> 

2800  GOSUB  3958 

2810  Di  2*<  <U0-Ua>*"2+<  V0-Va>''2>yv.  3 
2820  K-K+l 

2830  IF  Di 2>*Di 1  THEN  2860 

2840  Xx<X,l)>X7 

2830  GOTO  2720 

2860  IF  K > 1  THEN  2890 

2870  X7*Xx  <  X , 1 )  -D 

2880  GOTO  2800 

2890  D»D'2 

2900  J*J+l 

2910  IF  J<  >15  THEN  2640 
2928  P-CLa»da0'Lai»dac  >~2 
2930  S»*<L**d*0-'<2*PI*L«>  >^2 

2940  De<X>-P-Ss*<Xx<X,2>~2-Xx<X, 1>~2>  !  Dc  IS  THE  DIELECTRIC  CONSTRNT 

2930  Epp<X>-Ss*2*Xx<X,2>*Xx<X, 1)  !  Epp  IS  THE  LOSS  FACTOR 

2960  T*-Epp<X>/Dc<X> 

2970  Tad<X)-Ta  i  Tad  IS  THE  LOSS  TANGENT 

2980  TadwD=a*<l-P)'Za 
2990  Cg“Al/'<2*Bl ) 

3000  Taws»Tadw*<  (Cgt-P^Dc  <X>  J/'CCg+P)  > 

3010  Tad<X>*Ta-Taws 
3020  Epp<X)*Dc  <X)*Tad<X> 

3030  HEXT  X 

3040  IF  El«2  THEN  3070 

3090  IF  E>0  THEN  2350 

3060  PRINTER  IS  0 

3070  PRINT 

3080  PRINT 

3090  PRINT  "DIELECTRIC  CONSTANT";"  LOSS  FACTOR". "  LOSS  TANGENT" 

3100  PRINT 

3110  FOR  X»i*E  TO  3+E 

3120  IMAGE  4X.DDDD.DDD, 18X.MD.DDE, 10X.MD.DDE 
3130  PRINT  USING  3 l 20; Dc < X ) , Epp < X > , Tad < X > 

3140  NEXT  X 
3130  PRINTER  IS  16 

3160  INPUT  “DO  YOU  UANT  SOME  OTHER  POSSIBLE  VALUES  CY  OR  N]",Y* 

3170  IF  Y*-“Y*  THEN  2040 

3130  OUTPUT  9; "R*qu*st  tii»*“  !  FROM  THE  98033A  REAL  TIME  CLOCK 

3190  ENTER  9; Month, Day , Hour , Mi  nut*, Second 
3200  PRINTER  IS  0 

3210  PRINT 

3220  GOSUB  3910  !  SUBROUTINE  TO  HELP  PRINT  OUT  THE  TIME  UNDERSTANDABLY 

3230  PRINT 

3240  IMAGE  "Th*  n»«  »t  th*  *nd  of  tht  t**t  is:  " ,  22,  A,  22,  A,  ZZ 
3230  PRINT  USING  3240 ; Hour , Co  1  * , M i nut  * , Co  1  * , S*c ond 

3260  PRINT 

3270  PRINTER  IS  16  1  RETURN  PRINTER  CONTROL  TO  CRT  SCREEN 

3280  LOCAL  719 

3290  CLEAR  7 

3300  END  '  END  OF  MAIN  PROGRAM 


A-6 
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3310  !  - 

3320  !  THIS  SUBROUTINE  ZEROS  IN  ON  THE  VALUE  OF  dBm  THAT  CORRESPONDS 

3330  !  TO  A  VOLTAGE  READING  GREATER  THAN  100  mV 

3340  OUTPUT  7, 19; “PL" , Dbm, “DM“  !  SET  POWER  ON  HP-3330A  SWEEP  OSCILLATOR 

3350  TRIGGER  7,23  1  ASK  HP-3438A  METER  TO  OUTPUT  ITS  READING 

33S0  ENTER  7, 23; Volt*  !  READ  VOLTS  FROM  METER 

3370  LOCAL  7,23 

3380  IF  Vo  1 1  *  > . 08  THEN  GOTO  3430 

3390  Dbm»Db«+D*t«pl 

3400  OUTPUT  7, 19? *PL“,Dbm, “DM" 

3410  TRIGCER  7,23 
3420  ENTER  7, 23; Volt* 

3430  LOCAL  7,23 

3440  GOTO  3380 

3430  IF  Vo  1 1 * >. 095  THEN  GOTO  3320 

3480  Dbm»Dbm+D*t*p2 

3470  OUTPUT  7, 19; "PL", Dbm, “DM* 

3480  TRIGGER  7,23 
3490  ENTER  7, 23; Volt* 

3300  LOCAL  7,23 
3310  GOTO  3450 

3320  IF  Vol t*>*.  1.  THEN  GOTO  3390 

3330  Db**DbM*D*tep3 

3340  OUTPUT  7,  19;  "PL",  Dbm,  “DM" 

3530  TRIGGER  7,23 
3580  ENTER  7, 23; Volts 
3570  LOCAL  7,23 
3380  GOTO  3520 
3590  RETURN 

3600  !  SUBROUTINE  TO  MOVE  THE  PROBE  BY  COMMANDING  THE  KLINGER 

3610  EOL  7; CHR*< 13)  !  SET-UP  NEW  END  OF  LINE  STRING 

3620  OUTPUT  7,7  USING  “#,K,L*;"R  233*  ! 

3630  OUTPUT  7,7  USING  "#,K,L";*S  2"  .  1  R,$,F  SET  MOVE  SPEED  OF  PROBE 

3640  OUTPUT  7,7  USING  “*,K,L“;“F  20“  ! 

3630  OUTPUT  7,7  USING  " t , K, L" ; NsttpS  !  MOVE  THIS  NUMBER  OF  STEPS 

3660  OUTPUT  7,7  USING  “*,K,L";Pm*  1  DIRECTION  OF  MOVE 

3670  OUTPUT  7,7  USING  "#,K,L“;"G"  !  GO... EXECUTE  MOVE 

3680  LOCAL  7,7 
3690  RETURN 

3700  !  THIS  SUBROUTINE  DEFINES  THE  STEP  SIZE  MOVE  COMMAND  FOR  THE  KLINGER 

3710  N*«VAL*<St«p) 

3720  Q*«"N  " 

3739  Nst«p*«QSlN* 

3740  RETURN 

3730  !  THIS  ROUTINE  FINDS  THE  3  dBm  POSITION,  STARTING  AT  A  MINIMUM 

3760  I-St«p 

3770  GOSUB  3600  !  MOVE  PROBE 

3730  WAIT  300 

3790  TRIGGER  7,23 

3800  ENTER  7, 23; Volt* 

3810  IF  <Volt»«.l)  OR  <Volt*>.l)  THEN  GOTO  3840 

3820  I»I*St«p 

3830  GOTO  3770 

3840  LOCAL  7,23 

3838  RETURN 

3860  1  SUBROUTINE  TO  ASK  THE  KLINGER  TO  OUTPUT  THE  POSITION  OF  THE  PROBE 
3870  EOL  7 ; CHR*  < 13) 

3880  OUTPUT  7,7  USING  " » , K , L ’ 5 “ V  0"  '  OUTPUT  POSITION  FROM  KLINGER 

3890  ENTER  7,29;Po*  !  READ  POSITION 

3900  RETURN 

3910  !  THIS  SUROUTINS  HELPS  US  PRINT  OUT  THE  DATE  t,  TIME  IN  GOOD  FORM 

3928  Sl**h*»"/“ 

3930  Col *■“ 1  * 

3940  RETURN 

3930  *  THIS  IS  A  SUBROUTINE  FOR  THE  DIELECTRIC  CONSTANT  PART  OF  THE  PROGRAM 

3960  S»X6«2 
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3970 

3980 

3990 

4000 

4010 

4020 

4030 

4040 

40S0 

4060 

4070 

4080 

4090 

4100 

4110 

4120 

4130 

4140 

4130 

4160 

4170 

4100 

4190 

4200 

4210 

4220 

4230 

4240 

4230 

4260 

4270 

4^80 

4290 

4300 

4310 

4320 

4330 

4340 

4330 

4360 

4370 

4380 

4390 

4400 

4410 

4420 


T  »X7  *2 

Vl«X6*FNSi nh<S>*X7*SIN<T> 

V2*X6*SlN<T)-X7*FNSi nh<S> 

Bot  •  <  X6~  2fX7~2  >  *<  FNCos*<  S  > +C0S 

U«-Vl/Bot 

V*-V2/Bot 

RETURN 

!  SUBROUTINE  TO  FIND  THE  NEXT  ESTIMATE  OF  MINIMUM  POSITION 
N»1 

COSUB  3320  !  ZERO  IN  ON  THE  dBm  LEVEL 

OUTPUT  7,19; “OPPL"  !  ASK  SWEEP  OSCILLATOR  ITS  POWER  LEVEL 

ENTER  7,19;Db»<N>  !  READ  POWER  LEVEL  IN  dB» 

GOSUB  3860  !  ASK  THE  KLINGER  THE  POSITION  OF  THE  PROBE 

Pos<N)«Pos 
Vol ts<N)"Vo1 tt 
Db«"Db«-Dbl *ss 


Db«"Db«-Dbl*ss  !  SET  POWER  LOWER  BEFORE  NEXT  MOVE 

!  PRINT  LINES  ARE  COMMENTED  OUT  BUT  ARE  AVAILABLE  FOR  USE  IF  NEEDED 
!  PRINT  N;*  Position  »" ;  Pos<N) ;  “  <u»>  “ ,  “dB*  ■  Db*<N> ,  “Volts  *";VoU*(N> 
OUTPUT  7, 19; "PL", Db«, “DM" 

IF  N>-2  THEN  GOTO  4210 

N»N+1 

GOSUB  3600  .  !  MOVE  PROBE 

WAIT  1000  !  WAIT  1  SECOND  FOR  PROBE  TO  SETTLE 

GOTO  4060 

L«N 

N 1 *N—  1 

N-N+l 

Db»"Dbn<L>-Db 1 •*% 

OUTPUT  7, 19; "PL"  , Db«, “DM" 

GOSUB  3600  !  MOVE  PROBE 

WAIT  1000.  !  WAIT  1  SECOND  FOR  PROBE  TO  SETTLE 

GOSUB  3320  !  ZERO  IN  ON  dB»  LEVEL 

OUTPUT  7,19; "OPPL“ 

ENTER  7,19; Db*(N) 

GOSUB  3860  !  ASK  THE  KLINGER  THE  PROBE'S  NEW  POSITION 

Pos<N)"Pos 

Volts<N>"Volts 

!  PRINT  n;"  Position  -" ; Pos<N> ;  " <u«o " ,  *dBi«  ■  ■ ;  Db»CN> , "Vol  ts  -*;Vo1ts 
IF  DbM<N><-23. 0  THEN  GOTO  4210 

IF  <Obn<N)<>Dbm<L>)  AND  (Db«<L) >«Db»<Nl > >  THEN  GOTO  4380 
GOTO  4210 

p»S«"-« 

GOSUB  3600  !  MOVE  PROBE  BACK  TO  LOCATION  OF  MINIMUM  ESTIMATE 

Dbm«Db*<L> 

!  PRINT 
RETURN 


!  MOVE  PROBE 

!  WAIT  1  SECOND  FOR  PROBE  TO  SETTLE 
!  ZERO  IN  ON  <Stm  LEVEL 


!  ASK  THE  KLINGER  THE  PROBE'S  NEW  POSITION 


’ ; Db»(N) , -Vol ts  ""{Volts 
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APPENDIX  3 


DIELECTRIC  CONSTANT  CALCULATION  PROGRAM 


B-l 
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le  !  JULY,  1983  R.  WATKINS 

1  1  ! 

20  !  THIS  PROCRAM  IS  BASED  ON  THE  NOV.  1972  REPORT  AR3-NC-4  FROM  THE 

30  !  AGRICULTURAL  RESEARCH  SERVICE.  IT  WAS  REVISED  FROM  THE  PREVIOUS 

40  !  VERSION  WRITTEN  IN  1980  AND  NOW  IS  ACCURATE  FOR  ALL  CASES, 

j  30  !  IT  WILL  FIND  DIELECTRIC  CONSTANT,  LOSS  TANGENT, 

60  !  AND  LOSS  FACTOR  AND  WILL  PRINT  THEM  ON  THE  PRINTER  OR  THE  SCREEN. 

[70  •  AVG.  RUNNING  TIME:  1-3  MIN. 

80  !  DETAILED  INSTRUCTIONS  ARE  OBTAINED  BY  TYPING  'EDIT  1680'. 

90  !  MAKE  SURE  THAT  WAVECUIDE  DIMENSIONS  (LINES  1640-1660'  ARE  CORRECT 

100  !  BEFORE  RUNNING  PROGRAM. 

110  DEF  FNSinh<X)-.3+<EXP<X>-£XP<-X>> 

120  DEF  FNCosH <X ) ■ . 5*<EXP<X >  +EXP < -X )  ) 

!  130  DIM  Xx<9,2>,Dc<9>,Epp<9>,T*d<9),F<3,30>, Samp  1  SC 60 3 

•  140  GOSUB  1640 

‘  130  LINPUT  “SAMPLE  DESIGNATION* , S*apl * 

1160  INPUT  "FR£QUENCY<GHz)“,F 

170  INPUT  “LENGTH  OF  SAMPLE  <  i»»  >*,  L* 

180  INPUT  “2  OdB  POINTS  WITH  SAMPLE <»«>", Lzs , Rz» 

I  190  INPUT  "SAMPLE  Db“,Sd 

|  200  INPUT  “2  OdB  POINTS  WITHOUT  SAMPLE <mm>< IF  NONE,  INPUT  0,0",LZ4,Pz* 

|  210  INPUT  “dB  LEVEL  WITHOUT  SAMPLE* , Db 

,  220  INPUT  “ENTER  SHIFT”, W1 

>  230  M-X-0 

!  240  L**o«299.6966'F 

230  Luc>2«J 

260  LMgaLuO''(l*(Luo/L4tc)>'2)A.S 
j  270  IF  Lz*<>0  THEN  300 

Y  280  Lz**L*»g^2 

k’  290  Rza-LMg/2 

[  300  Sc*l+<W*L*»g>^2'<8*PI*B1 *A1 -3) 

L  310  2**. 5*<Lz*+Rz*) 

I  320  Zl«Za 

P  330  <3-0 

J.  340  Q-Q+l 

t  330  IF  ABS<Q*<L*»g/'2)-<Z*+L*«g/2)  ><L**g^4  THEN  370 

S360  GOTO  340 

370  Zufl*(Lug/2) 

380  Zs-. 3*<Lz**Rzs> 

390  Dz*«ABS<Lz*-Rz*)/Sc 
400  Dzs-ABS<Lzs-Rz*)'Se 
416  St-<Zl-Zs>'Se 

I  420  IF  Dz»-0  THEN  460 

J  430  Apg*. 2302383*Db 

,  440  Po-EXP(Arg) 

j  430  Sur-<Po-C0S<PI*Dz*/L*»g)-2)-. 3^SIN<PI*DZ4/L*«g) 

I  460  Arg-.2302383*Sd 

1  470  Po-EXP< Arg) 

P  480  IF  Dza-0  THEN  320 

(  490  Arg«< < Po- l) ^ < Sur-2- l ) )-. 3 

I'  300  Dz**L*«g*ASN<Arg)/PI 

310  R-<Z*-L*-3t  "Z* 

320  Ppr-P  I  •*<  Dz»-R*Dza;  ^L**g 

330  Sur«<Po-COS(Ppr;A'C^'‘.  5^SIN<Ppr> 

340  Phi-l/Swr 
330  N-0 

360  IF  - <  L**St  ) *.  3*L*n»g*N>0  THEN  390 
i,  370  N-N+l 

f  380  GOTO  360 

r  390  20--<3t*L*.>*.  3*L*mg»N 

L  o>  .  U*2*P  I  *Z0''L*i»g 

•f  10  •i*TAN<U^»<'Pbi  a2-1  ) 

S»Phi  *(  1*TAN<U)a2) 

t'O  C*',A'S2*B'2>A.  3*<L*(i>'3  i  *■  P  I  ■*  L  *  •*  '  1  *  <  P  h  i  *  T  A  N  (  U  J  |AI)  >  ■ 

^>•*0  C»C*<  A/ ABS  <  A  ) 

:3<?  IF  C<1  THEN  K-2 


660  FOR  Y«K*-1*W1*2  TO  6-»K*Ul*2  STEP  2 

670  Mt«Y*PU2 

630  H2«<  Y-2>*PI'2 

690  IF  M2<8  THEN  n;»0 

700  X2-<(1t-(12)  -2*fl2 

710  Fu»TflN(X25  'X2-C 

720  IF  Fu>0  THEN  Ml«X2 

730  IF  Fu< *0  THEN  PI2-X2 

740  IF  II  1-H2<.  000000001  THEN  760 

730  GOTO  700 

760  H«M*1 

770  Xx<(1,  n»X2 

780  ftd»Xx<H,  1  >*l_4pg*PHi  ■*<  1  ♦TflN< U ) "2  ) 

790  B«n»2*PI*L**<l*TftN<;<x<M, 1  > >-2-TRN<XxCM, l>)^Xx(M, l,, 

800  Xx<M,2)*Ad/'D*n 

810  NEXT  Y 

820  cf*L**g'<2*PI*L«> 

830  l*<Phi  *TflN<U)  )^2) 

840  U0»A»D»n 
830  V0«-B*D«n 

860  FOR  X-l  TO  3 
870  J-0 

880  X6>Xx(X,2) 

890  X7-Xx<X,l) 

900  O-.01 
910  GOSUB  1360 

920  Di l*( <U0-U*)^2*<VO-V*>^2)^. 3 
930  IF  Dil<. 00001  THEN  1220 
940  K»0 

930  X6«Xx<X,2)+D 

960  X7»Xx<X,l> 

970  GOSUB  1360 

980  Di2»<<U0-U4>-2*<V0-V4>-2>«.3 

990  K-K+l 

1000  IF  D12>-Dil  THEN  1040 

1010  Xx<X, 2 >*X6 

1020  Oi 1-Di2 

1030  GOTO  930 

1040  IF  K > 1  THEN  1070 

1030  X6«Xx<X,2)-D 

1060  GOTO  970 

1070  K«0 

1080  X7»Xx<X, 1 )+D 
1090  X6«Xx<X,2) 

1100  COSUB  1360 

1110  DI  2«<  <UO-U*>  ~2*<  V0-V*;-'2)-<  3 

1120  K-K*l 

1130  IF  Di 2>aDi l  THEN  1160 

1140  Xx  <  X , 1 )»X7 

1130  GOTO  1020 

1160  IF  K > 1  THEN  1190 

1 1 70  X7*Xx ( X , 1 )-D 

1130  GOTO  1100 

1190  D-D'2 

1200  J«J*l 

1210  IF  ;<>U  THEN  940 

1220  P»< l*»0''L**c  >~2 

1230  S*«<L*»o'<2*PI*L* i '"2 

1240  De <X)*P-S»*<xx<x. 2;~2-Xx<x, 1 j-2) 

1230  Epp<x>«s**2*xx<:<,  2>-xx<::,  i  - 
1260  T*«Eppex  •✓De ''X‘> 

1270  T»g<X)»T4 

1280  T4dw»Dz4*<  l-P>  -  Zh, 

1290  Cg»«l  >•  C  2* B  T  ) 

1388  T*«»»T*du*<  (Cg-»P  'Dc  Cw  )-'<Cg*P>  > 

1310  T*d<;<)«T*-T*w* 
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1320 

1330 

1340 

1330 

1360 

1370 

1380 

1390 

1400 

1410 

1420 

1430 

1440 

1441 
1430 
1460 
1470 
1480 
1490 
1300 
1310 
1320 
1330 
1540 
1330 
1560 
1570 
1580 
1390 
1600 
1610 
1620 
1630 
1640 
1630 
1660 
1670 
1680 
1690 
1700 
1710 

1720 

1721 

1721 

1722 
1730 
1740 
1730 
1760 
1770 
1780 
1790 
1300 
1310 
1820 

1330 

1331 
1832 
1340 
1830 
I860 
1870 
1880 
1881 
1330 


cppvX-i*DcvX>*T*d(X> 

NEXT  X 

IF  El-2  THEN  1360 
IF  E>0  THEN  658 

INPUT  “PRINTER  OR  SCREENS  OF  S>“,Y* 
IF  Y*-“P*  THEN  PRINTER  IS  0 
PRINT 

PRINT  “SAMPLE  DESIGNATION; 

PRINT 

'FREQUENCY  IN  GIGAHERTZ 


' ; 3  sap 1  * 


:  f 


•sample  length  in  millimeters:  ";l« 


•SHIFT«“;wi 


LOSS  FACTOR  * 


PRINT 
PRINT 
PRINT 
PRINT 
PRINT 
PRINT 
PRINT 

PRINT  “DIELECTRIC  CONSTANT* : ' 

FOR  X-l*£  TO  3*E 

PRINT  Dc <X) , Epp<X) , TsdCX) 

NEXT  X 
PRINT 

PRINTER  IS  16  • 

INPUT  “ANOTHER  VALUE  OF  SHIFT < Y  OR  N>",YS 

IF  Y*-“Y“  THEN  220 

STOP 

S«X6*2 

T -X7*2 

V 1 -X6*FNS i nh < S > ♦X7*S I N < T ) 
V2«X6*SIN<T)-X7*FNSinh<S> 
Soi-<X6-2+X7-2>*<FNCo*h<S>+C0S<T) ) 

Ua-V 1 /Bot 


LOSS  TANCENT" 


V4*V2/'Bot 
RETURN 
A1 -22.8680 
B1 -10.  16 


*  UIDTH  OF  WAVEGUIDE 
HEIGHT  OF  WAVEGUIDE 


IN  MILLIMETERS 
IN  MILLIMETERS 


W-3.2 

RETURN 


!  SLOT  WIDTH  IN  MILLIMETERS  OF  THE  WAVEGUIDE  SLOTTED  LINE 


INSTRUCTIONS  FOR  USE  OF  DIELECTRIC  CONSTANT  PROGRAM 


R.  UAtkins 


1.  If  you  don't  need  to  enter  wAveguide  cross  sectionAl 
dimensions,  go  to  2. 
dimensions,  go  to  2. 


a.  Push  EDIT  1640  if  you  need  to  ehsnge  wAveguide  dimensions 

b.  Use  cursor  controls  to  position  cursor  on 
the  line  thst  needs  to  be  chsnged. 

c.  Use  cursor  controls  to  position  cursor  on 
the  number  thst  needs  to  be  chsnged 

d.  Enter  chsnges. 

e.  When  line  is  correct,  push  STORE  And  go  on  to  next  line. 

f.  When  sll  is  correct,  push  PAUSE  to  exit  the  edit  mode. 


Push  RUN 


Enter  the  S  Ample  Designs*  ton  •.  desc  r  t  pt  i  on  of  ssmoie.',  then  CONTINUE 
Enter  frequency  in  gigshert:,  then  CONTINUE. 

Enter  ssmple  length  in  mm,  then  CONTINUE. 


Enter  tuo  positions  Along  the  slotted  line  in  mm  where  the 
power  i s  up  s  given  number  of  dB  f rom  the  minimum  dB  vsl ue 
with  the  ssmple  in.  These  points  should  be  messured  relstiwe  to 


8-4 
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the  »hort.  Separate  vtlu*l  Dy  4  coaat,  th*n  CONTINUE. 

7.  Enter  the  dB  level  of  th«  xuurtarmi  in  iuc  a. (usually  2) 

*nd  then  CONTINUE 

8.  3m«  as  step  8.  but  with  th*  waveguide  empty.  Then  CONTINUE. 

9.  dB  level  of  *«*sur«Mnt  in  step  3.  Then  CONTINUE. 

10.  Enter  shift.  This  is  4  positive  whole  number  th4t  determines 

whst  3  possible  values  of  dielectric  constant  will  be  examined 
Then  CONTINUE.  If  you  have  no  idea  what  value  to  use  try  9, 
then  1 


When  asked,  enter  P  (printer)  or  S  (screen)  to  determine 


where  data  will  be  displayed 


APPENOIX  C 


SAMPLE  COMPUTER  OUTPUT 


Date:  04/10/85 

Time:  11:55:27  At  start  of  test 

Relative  Humidity  *  50  percent  Temperature  *  75  Deg.  F 

Frequency  selected  for  this  test  is:  8.400  (GHz) 

Waveguide  Wavelength  *  57.055  (mm) 

THIS  IS  A  TEST  OF  THIS  PROGRAM  TO  SEE  IF  IT  WILL  GIVE  DC-1  FOR  AIR 
Length  of  Sample  is:  25.4  (mm) 


DIELECTRIC  CONSTANT 


LOSS  FACTOR 


LOSS  TANGENT 


1.000 

2.207 

4.297 


-1.72E-04 

-2.28E-03 

-3.22E-03 


1.72E-04 

1.03E-03 

7.49E-04 


The  time  at  the  end  of  the  test  is:  12:05:58 


Oate:  04/10/85 

Time:  14:49:38  At  start  of  test 

Relative  Humidity  «  50  percent  Temperature  -  80  Deg.  F 

Frequency  selected  for  this  test  is:  3.400  (GHz) 

Waveguide  Wavelength  -  57.055  (mm) 

THIS  IS  A  TEST  OF  THIS  PROGRAM  TO  SEE  IF  IT  WILL  GIVE  DC-1  FOR  AIR 
Length  of  Sample  is:  12.7  (mm) 
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DIELECTRIC  CONSTANT 


LOSS  FACTOR 


LOSS  TANGENT 


1.000 

4.950 

12.843 


2.06E-04 

-5.67E-03 

-1.70E-02 


2.06E-04 

-1.15E-03 

-1.32E-03 


The  time  at  the  end  of  the  test  is:  15:01:11 
Date:  04/01/85 

Time:  15:4  At  start  of  program 

Relative  Humidity  *  56  percent 

Frequency  selected  for  this  test  is:  8.3  (GHz) 
Waveguide  Wavelength  *  58.864  (mm) 

WHITE  CERAMIC  FOAM 

Length  of  Sample  is:  9/32  (mm) 


Temperature  a  82  Deg.  F 


DIELECTRIC  CONSTANT 

1.507 

9.012 

24.02 


LOSS  FACTOR 

2.11E-03 

-8.32E-03 

-2.89E-02 


LOSS  TANGENT 

1.40E-03 

-9.23E-04 

-1.20E-03 


The  time  at  the  end  of  the  test  is:  15:25 


Date:  04/19/85 

Time:  14:27:40  At  start  of  test 

Relative  Humidity  *  43  percent 

Frequency  selected  for  this  test  is: 
Waveguide  Wavelength  *  57.055  (mm) 

POLYSULFONE  POLYMER  SAMPLE 

Length  of  Sample  is:  12.903  (mm) 


Temperature  *  79  Deg.  F 


8.400  (GHz) 


DIELECTRIC  CONSTANT 

2.971 

9.856 

20.942 


LOSS  FACTOR 

2.97E-02 

8.44E-02 

1.27E-01 


LOSS  TANGENT 

1.00E-02 

8.57E-03 

6.05E-03 


The  time  at  the  end  of  the  test  is:  14:44:46 
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Department  of  the  Navy 
Washington,  DC  20361 


Commander 

Naval  Sea  Systems  Command 
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Washington,  DC  20362 

Oi rector 

Strategic  Systems  Program 
Office  (PH-1) 

Attn:  SP2731  (J.  Culver)  I 

SP273  (E.  Throckmorton)  1 

Department  of  the  Navy 
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1 
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1 
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1 
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1 
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1 
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1 
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1 
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1 
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1 
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1 
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10 
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(D. 

Price) 

1 

R13 

(S. 
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